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Abstract 

We investigate the quantitative impact that data from the second oscillation maximum has on the perfor- 
mance of wide band beam neutrino oscillation experiments. We present results for the physics sensitivities to 
standard three flavor oscillation, as well as results for the sensitivity to non-standard interactions. The quan- 
titative study is performed using an experimental setup similar to the Fermilab to DUSEL Long Baseline 
Neutrino Experiment (LBNE). We find that, with the single exception of sensitivity to the mass hierarchy, 
the second maximum plays only a marginal role due to the experimental difficulties to obtain a statistically 
significant and sufficiently background-free event sample at low energies. This conclusion is valid for both 
water Cerenkov and liquid argon detectors. Moreover, we confirm that non-standard neutrino interactions 
are very hard to distinguish experimentally from standard three-flavor effects and can lead to a considerable 
loss of sensitivity to #13, the mass hierarchy and CP violation. 
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I. INTRODUCTION 



Neutrino physics has seen a spectacular transition from a collection of anomalies to a field of 
precision study with a firmly established theoretical underpinning. All neutrino flavor transition 
data, with the exception of LSND pQ and MiniBooNE [2j, can be described by oscillation of three 
active neutrino flavors, see e.g. the reviews in [3JU]. Throughout this paper we assume that LSND 
and MiniBooNE have explanations which do not affect our results, i.e. they are not due to neutrino 
oscillation. 

Within the three flavor oscillation framework, experiments have determined the values of the 
mass squared differences and the associated large mixing angles with a precision at the level of a 
few percent [3J. Currently unknown are the size of #13, the value of the CP phase 5qp and the sign 
of the atmospheric mass splitting Am^, as well whether 623 is exactly vr/4, and if not, whether it 
is larger or smaller than 7r/4. The ultimate goal is to determine the neutrino mixing matrix with at 
least the same level of precision and redundancy as the CKM matrix in the quark sector. The size 
of #13 plays a particularly crucial role, since this quantity will set the scale for the effort necessary 
to answer the open questions. The need to determine #13 has spurred a number of reactor neutrino 
experiments [5] using disappearance of v e : Double Chooz [6j, RENO |7J and Daya Bay [8]. Their 
discovery reach at 3d confidence level will go down to approximately sin 2 2^13 = 10 -2 [HI [10]. At 
the same time the next generation of long baseline experiments looking for Va —> v e is underway 
with T2K and NOz^A |12j . Neither T2K nor NOi^A can provide information on <5cp beyond 
a mere indication and even that is only possible in combination with the data from Daya Bay [9]. 
The discovery of the mass hierarchy, if discovery is defined at the usual 5 a confidence level, will 
not be possible; even a 3 a evidence is very unlikely [9] . Therefore, it has been widely recognized 
that long baseline experiments with physics capabilities far beyond NO^A and T2K are necessary, 
for a review of the possibilities, see reference |13j . 

Here, we would like to focus on the superbeam concept, or more specifically on what is called 
a wide band beam (WBB). In a superbeam, muon neutrinos are produced by the decay of pions, 
where the pions have been produced by proton irradiation of a solid target. All neutrino beams 
relevant in the context of this study use a magnetic horn to focus and sign-select the pions. In a wide 
band beam, the detector is on-axis and thus receives a wide (sicl) energy spectrum of neutrinos. 
The wide band beam concept makes maximal use of the available pions and thus provides higher 
event rates compared to a narrow band or off-axis beam. The price to pay for the wide spectrum is 
that the detector needs to have a very good energy resolution, and the existence of a high energy 
tail in the beam will lead to feed down of neutral current background events. Thus, a wide band 
beam imposes unique demands on the detector technology. Apart from allowing for more events, 
the wide beam spectrum allows to study a range of L/E values within one experiment, and possibly 
even to observe more than only one oscillation maximum^] On the level of oscillation probabilities, 
the ability to observe two or more cycles of the oscillation obviously allows to distinguish between 
otherwise degenerate solutions. Therefore, the observation of the second oscillation maximum is 
considered to play an important role in wide band beam experiments. The purpose of the present 
paper is to study in detail and in a quantitative manner whether the assertion of the role of the 
second oscillation maximum based on probabilities remains valid in a full numerical sensitivity 
calculation. We also include the case of non-standard interactions, where one expects similar 
benefits from the presence of the second oscillation maximum. 

In order to perform a full numerical sensitivity calculation we need to specify the experimental 
parameters in great detail and therefore have to constrain the numerical analysis to a specific setup, 
which we model to resemble the Fermilab to DUSEL Long Baseline Neutrino Experiment (LBNE). 

1 We will use the term oscillation maximum also for disappearance channels, where actually a minimum in the 
survival probability is observed. 



2 



However, whenever the specifics of the chosen experimental setup obscure the underlying physics, 
we will show results for sensible variations around our chosen setup. In section [TT] we discuss the 
theoretical framework with respect to standard and non-standard oscillations. In section III we 



spell out the details of the experimental setup and describe our analysis techniques. Section IV 



will contain our results on both three flavor oscillation and non-standard interactions and finally 
in section [V] we will summarize our findings and present our conclusion. In appendix [A] we show 
supplementary results on variations of the total exposure and baseline. 



II. FRAMEWORK 



A. Three flavor oscillation 



In this paper our main concern is the measurement of the transition probabilities P{yu v e) 
and P(Pfj, — > v e ). Since the baseline considered is longer than 1, 000 km, matter effects will play an 
important role. While the underlying Hamiltonian describing oscillations in the presence of a matter 
potential is quite simple, the resulting expressions for the exact oscillation probabilities are not. 
Therefore, a plethora of approximations have been devised, for an overview see reference |14j . Even 
these approximate solutions have a very rich structure; in particular, for fixed energy, any given 
value of the oscillation probability can typically be realized by different combinations of oscillation 
parameters. The possible degenerate solutions can be classified into the intrinsic ambiguity |15j , 
the sign of Am^ ambiguity [16] and the octant ambiguity [17] . Combined, these three ambiguities 
give rise to what has become known as the eightfold degeneracy [18] . A recent comprehensive 
analytical discussion of the eightfold degeneracy can be found in reference [19]. In the context of 
long baseline experiments the most worrisome degeneracy stems from a combination of the intrinsic 
and sign ambiguity which can lead to a phenomenon called 7r-transit [20J, in which a CP violating 
true solution is mapped into a CP conserving fake solution. A large number of possible remedies 
has been proposed in the literature, here we focus on the proposal to use not only neutrino and 
anti-neutrino events from the 1 st oscillation maximum, but also from the 2 nd oscillation maximum. 

In figure [l] we illustrate how the use of the 2 nd oscillation maximum can alleviate the sign 
degeneracy. At a baseline of 1,300 km, the first oscillation maximum occurs with Am^ = 2.4 x 
10 _3 eV 2 for E v = 2.5 GeV and the second oscillation maximum is at E v = 0.84 GeV. The first 
zero of the oscillation term occurs at E v = 1.25 GeV, and we will use this energy to separate events 
from the two oscillation maxima, assigning all events below to the second maximum and all above 
to the first. In figure [T] we show so called bi-rate plots. In a bi-rate plot #13 is kept fixed and the 
coordinates are the total number of events in the neutrino channel and in the anti-neutrino channel, 
respectively. For each possible choice of the CP phase one obtains a point in this kind of diagram, 
and as the CP phase is continuously varied from — ir to +ir the points trace out a banana-shaped 
curve (or, on a linear scale, an ellipse) [21]. This is similar to a bi-probability plot [16] . but avoids 
the problem of choosing a neutrino energy for plotting and thus allows a closer approximation of 
the experimental realities. In figure [TJ we moreover separated the event sample into events below 
1.25 GeV (green/light gray curves) and events above that energy (blue/dark gray curves). The solid 
lines are for normal hierarchy and the dashed ones for inverted. The solid disks indicate the event 
rates for 5cp = — 7r/2, whereas the open circles indicate the event rates for 5qp = +vr/2. Focusing 
on the 2 nd maximum (green/light gray lines), we see that the "bananas" for both hierarchies are 
very similar and occupy essentially the same area in the event rate plane. Moreover, the event 
rates at the two maximally CP violating values of <5cp = ±7r/2 do not change when going from 
one hierarchy to the other. For the 1 st maximum (blue/dark gray lines), the two "bananas" are 
very different and the event rates for the same value of <5cp change greatly when switching the 
hierarchy. Therefore, given enough statistics, we expect the measurement in the 2 nd maximum to 
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FIG. 1. Bi-ratc plot for a typical wide band beam for the 1 st (blue/dark gray) and 2 nd (green/light 
gray) oscillation maximum for normal (solid) and inverted (dashed) hierarchy. Along each of the banana- 
shaped curves, the CP phase i5cp varies from —it to tt, where <5cp — —tt/2 is denoted by a solid disk and 
5qp = +7r/2 by an open circle. We have assumed sin 2 26*i3 = 0.025, and the other oscillation parameters 
are chosen according to eq. 

provide a clean value for <5cp, but no information on the hierarchy. The measurement in the 1 st 
maximum, on the other hand, should yield strong evidence for the mass hierarchy, but may suffer 
from degeneracies for the determination of 5qp ■ The combination of the two maxima should result 
in a clear and unambiguous determination of both the mass hierarchy and the CP phase. 

B. Non-standard Interactions 

At energies of a few GeV, relevant to accelerator neutrino oscillation experiments, the effects of 
new physics, which is expected at or above the electroweak scale, can be parametrized in terms of 
an effective theory. Some types of low-scale new physics can also be parametrized that way |22[ [23] . 
A well known example for the use of effective theory is the Fermi theory of nuclear beta decay. 
In this paper, we will use such non-standard neutrino interactions (NSI) as a benchmark scenario 
for deviations from the standard three-flavor oscillation framework, but we should keep in mind 
that new physics in the neutrino sector can also have different manifestations; examples are CPT 
violation or mixing between active and sterile neutrinos. Typical operators inducing non-standard 
neutrino interactions (NSI) are 

£cc => -2V2G F e^ JJ ' [v a ~i p P L lp\ [h P P L f] + h.c. (1) 

(charged current NSI) and 

£nc => -2V2G F e^ J [va-fPfp] [/VV] + h.c. . (2) 

(neutral current NSI). Here, Gf = V^g 2 /8Myy is the Fermi constant, Pl = (1 — 7 5 )/2, a and 
f3 are flavor indices of the neutrinos v and the charged leptons £, and the fermions / and /' are 
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CC f f' NC f 

the members of an arbitrary weak doublet. The parameters £ a p and £ a g give the relative 
magnitude of the NSI compared to standard weak interactions. For new physics around the TeV 
scale, we expect their absolute values to be of order 10 -3 -10 -2 . In the presence of new degrees of 
freedom below the electroweak scale, NSI could be larger and also expectations for the magnitude 
of NSI based on effective theory approaches [Ml [25] may be too conservative. Note that equations [T] 
and [2] include only V — A type interactions, but in principle, more general Lorentz structures are 
possible (see e.g. refs. [26, 27 \ for an overview). 

For phenomenological purposes, it is convenient to parametrize NSI in a slightly different way. 
Consider the u a — > i/g oscillation probability at baseline L, 



Pa(3 = \(vp\e- iHL Wa)\ Z , (3) 



with the Hamiltonian 



H = U\ Am 2 21 /2E + Vmsw , (4) 

V Am 2 31 /2Ej 

where U is the leptonic mixing matrix, E is the neutrino energy, and Vmsw is the 3x3 matrix 
describing matter effects. In the presence of CC NSI, the initial and final states get modified 
according to 

IO = K> + Yj £ U\ v p) and (^1 = M + Yj £ tp( u <*\ ' ( 5 ) 

f>=e,n,T a=e,fi,r 

respectively. The parameters e s aB and e^g, which are closely related to the parameters £ a g 
defined above [261 GZ] , describe non-standard admixtures to the neutrino states produced in associ- 
ation with a charged lepton of flavor a or detected in a process involving a charged lepton of flavor 
/3, respectively. Note that in e s a3 , the first index corresponds to the flavor of the charged lepton, 
and the second one to that of the neutrino, while in e^, the order is reversed. The matrices (l+e s ) 
and (1 + e ) need not be unitary, i.e. \v^) and are not required to form complete orthonor- 
mal sets of basis vectors in the Hilbert space. Instead of considering an oscillation probability 
P normalized to unity, it is therefore more useful to consider the apparent oscillation probability 
P( u a ~^ v< b)i defined as the number of neutrinos produced together with a charged lepton of flavor 
a and converting into a charged lepton of flavor (3 in the detector, divided by the same number in 
the absence of oscillations and non-standard interactions. The apparent oscillation probability is 
given by 



n<^4) = H\e- lHL \<)\ 2 



iHL\ ri 1 „s\ I 2 



;i + e % (e-*"%(i + 

[l + e d ) T e~^ L (l + e s ) T ] \ (6) 

J pa 

where H = t/diag(0, Am,2 1 /2£', Am 2 1 /2E) W + Vmsw- The modified matter potential is 

/l + £ m £ m m \ 
_ / ee e f- er \ 

Vmsw = V2G F N e e™ e™ , (7) 

\ p-m* pin* .m / 

with e™p being closely related to the £ a g from equation 2 As explained above, the magnitude of 
the e s ' d ' m parameters is expected to be at or below the 10 level for new physics at the TeV scale. 
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Proton energy 


pot per polarity 


Comment 


120 GeV 


22 x 10 20 


accelerator complex without Project X 


120 GeV 


72 x 10 20 


accelerator complex with Project X 



TABLE I. Neutrino beam configurations considered in our simulations. 



Model-independent experimental bounds on the s s,d,m parameters are typically of O(10~ 2 — 1) (3[ 
[28j 129] . In a specific model, however, the bounds may be much stronger because in most models 
neutrino NSI are accompanied by charged lepton flavor violation, which is strongly constrained 
by precision tests of the electroweak theory and by rare decay searches. Prom a model building 
point of view, it is therefore not easy to realize large non-standard neutrino interactions that can 
saturate the experimental bounds [2H [25] . 

Obviously, with any new experiment, we need to compare the expected bounds on new physics 
with the ones we already have. We use the bounds derived in [3[ l28l 29J as our benchmark. 
In particular, we use the 90% confidence level constraints |e™| < 4.2, |e™ | < 0.33, |e™| < 3.0 
\e™\ < 0.068, \e™ T \ < 0.063, |e™| < 0.2. Note that the relatively strong constraints on |e™ | and 
|e™ | have been derived from atmospheric neutrino data in a two- flavor framework; when three- 
flavor effects — in particular correlations between different types of NSI — are taken into account, 
these bounds may become somewhat weaker [30, 31J. 

III. METHODS 

A. Experimental setup 

To assess the sensitivity of a wide band neutrino beam to standard and non-standard oscillation 
physics, we have performed simulations using the GLoBES software |32|l33j. with an implementation 
of NSI developed in refs. [261 134| [35] . Our experiment description follows the LBNE proposal for a 
long-baseline neutrino beam from Fermilab to DUSEL, but our results will hold qualitatively also 
for other wide band beam experiments. 



1. Beam 

For the neutrino beam, we consider the options listed in table [I] We use the unit protons on 
target (pot) since it is the usual measure of integrated luminosity, C, for this kind of experiments. 
To compare this with other beams of different energy it is useful to convert this result to equivalent 
beam power, P, where we assume that the beam is on for 2 x 10 7 s per tropical year. 

P = 0.801 ( -^-) ( -== ^ kW. 

V GeV J V 10 20 pot per year J 

With this in mind, the luminosities given in table [I] correspond to about 6 years of running (3 
years in neutrino mode + 3 years in anti- neutrino mode) at a beam power of either ~ 700 kW 
or ~ 2,300kW. We have also studied the performance of a 60 GeV beam, which would have the 
advantage of lower backgrounds, at the expense of less statistics. Since we found only very minor 
performance differences between the 60 GeV and 120 GeV options, we restrict the discussion to the 
120 GeV beam in the following. Simulated spectra for all beam options have been kindly provided 
to us by the LBNE collaboration |36J. 
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2. Detectors 



We assume the far detector to be located at a baseline L = 1,300 km, corresponding to the 
distance from Fermilab to DUSEL, which translates into an energy E ~ 2.5 GeV for the 1 st oscil- 
lation maximum. In order to be sensitive to both oscillation maxima, and given the beam spectra, 
the detector needs to have good efficiency in the energy range from 0.5 — 4 GeV. The number of 
events in the 1 st oscillation maximum will be significantly larger than that in the 2 nd maximum. 
Currently, two detector technologies are considered in this context 

1. A water Cerenkov (WC) detector with a fiducial mass of 200 kt. As demonstrated by Super- 
Kamiokande, this type of detector allows for a clean separation of muon and electron quasi- 
elastic events. However, its application at GeV energies requires careful consideration of 
possible backgrounds from neutral current events giving rise to energetic neutral pions. A 
considerable amount of work went into studying this issue [371 EH] • Both studies agree quite 
well and we use the results from reference [37] . The GLoBES description of this WC detector 
is based on references |39|, 140] . This simulation includes energy-dependent efficiency tables, 
smearing matrices, and background estimates based on Monte-Carlo codes developed by 
the Super-Kamiokande collaboration |37j . We include both events from the v e appearance 
channel as well as the disappearance channel, for both neutrino and anti-neutrino running. 

2. A liquid argon (LAr) time projection chamber (TPC) with a fiducial mass of 34 kt. For the 
case of LAr, only much less detailed and accurate simulations are available since no large LAr 
detector has ever been operated. The LAr description is based on references [4H 142] . We 
include the v e and v e appearance channels, with backgrounds from the intrinsic f e /P e con- 
tamination of the beam, misidentified muons, and neutral current events. Our background 
estimate is conservative because the very high spatial resolution and the ability to detect 
very low energy particles in a LAr detector might allow for a much more efficient rejection of 
neutral current events. For the [p^ disappearance channel, the main backgrounds stem 
from neutral current events (we assume a rejection efficiency of 99.5%) and from the (y^) 
"wrong sign" contamination of the beam. Since the oscillation probabilities in the f« and 
Pp. disappearance channels are similar (except for the sub-leading contribution from matter 
effects), the latter background does not constitute a problem. 

For both detectors, the neutrino cross sections are based on [43^ 144] : they are computed for 
water and isoscalar targets, respectively. In an actual experiment great care needs to be taken to 
correctly model the cross sections, including nuclear effects. In our case, since we are using the 
same cross section to compute the data and perform the fit to that simulated data, any error due 
to the omission of nuclear effects will cancel. 

In our simulations we also include a near detector; for standard oscillation physics, its main 
effect is to reduce systematic uncertainties in the far detector, but for charged current NSI searches, 
it is valuable also as a standalone detector and its inclusion in the simulation is imperative. Since no 
specific technology has been chosen for the LBNE near detector(s) yet, we take a generic approach 
and assume the near detector to have identical properties (resolution, efficiencies, backgrounds, 
etc.) as the far detector, but a fiducial mass of only 1 kt. Furthermore, we assume the geometric 
acceptance of the near and far detectors to be the same, which greatly simplifies the calculation, but 
is very difficult to achieve in practice. In reality, the optimum choice of near detector technology and 
geometry may be very different for the WC and LAr cases, and thus also the effective systematic 
uncertainties may be quite different. 
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FIG. 2. Expected event rates in the v e appearance channel for sin 2 26*i3 = 0.05, (5cp = as a function of the 
reconstructed neutrino energy. The black histograms show the signal + background rates, while the filled 
red, blue, and beige histograms depict the backgrounds due to the intrinsic v e contamination of the beam, 
misidentified events, and neutral current events, respectively. Except for #13 and <5cp, the oscillation 
parameters are chosen according to eq. (pi). 



For illustration, we show in flg.[2]the expected event rates in the v e appearance channel for both 
detectors. It is clear that for small #13, backgrounds will be a limitation. In particular, neutral 
current events contaminate the second oscillation peak. 

B. Analysis 

To analyze the simulated data sets and to compute exclusion limits and allowed parameter regions, 
we use a x 2 analysis following ref. |20| . Our x 2 function has the form 



# of bins 2 (G\ _ C)( ^ 

\ 1 _ii 



X 2 =min[ E E xLt(^^^(©)>^+E7i]+E 

d=N,F S=V e ,V e ,Vn,i>p 3=1 i ^ 1 ' j ' 

(8) 



where Nffj and iV^ h SJ (0) are the observed and theoretically predicted event rates for detector d 

{d = N (near) or F (far)), event sample s (s = u e , v e , v^, u^), and bin j. The vector stands for 
the oscillation parameters, while a contains the systematical biases. The first term on the right 
hand side of equation [8] is the statistical contribution to x 2 > while the second term contains pull 
terms that disfavor values of the biases much larger than the associated systematic uncertainties 
erf. In a similar way, the last term of equation [8] is used to confine the oscillation parameters 
to within the region determined by other experiments, where, for each oscillation parameter Qj, 

0^ denotes the externally given best fit value and er® the la uncertainty on that value. In our 
simulations, we include such external prior terms only for the solar oscillation parameters 612 and 
Am|i to which the wide band beam is not sensitive. We assume the solar parameters to be known 
to within 5% at the la level, while for all other oscillation parameters, we set a® = 00. The default 
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wc 



LAr 



N/F correlated? 



Beam flux [%] 

Intrinsic background [%] 

Signal normalization for v e sample [%] 

Background normalization for v e sample [%] 

Signal normalization for sample [%] 

Background normalization for sample [%] 



oo 

OO 

0.7 
3.5 
0.7 
7.0 



oo 

00 

0.7 
7.0 
3.5 
7.0 



yes 
yes 
no 
no 
no 
no 



TABLE II. Systematic uncertainties assumed in our simulations. All systematic errors are assumed to 
be completely uncorrelated between the neutrino and anti- neutrino runs of the experiment. Note that 
uncertainties that are uncorrelated between the near (N) and far (F) detectors will add in quadrature when 
translated into an error on the measured oscillation probability. For example the near-far uncorrelated 
0.7% uncertainty in the number of v e signal events would translate into a 1% uncertainty on the measured 
oscillation probability. 



oscillation parameters used in this study are, in agreement with current fits [1], 



sin^ 12 = 0.32 , 0i 3 = , 

$23 = ~7 j ^CP = 7, 7T ■ 



+7.6 x 1(T 5 eV 2 , Amgi = +2.4 x 10~ 3 eV 



(9) 



and we use a conservative 5% uncertainty on the matter density. 

The systematic errors we include in our study are listed in Table |LT| for both detector technologies. 
We treat the normalization of the beam flux and of the background contributions to the v e and v e 
event samples as completely free parameters, i.e. we do not include pull terms for them. Moreover, 
we allow for uncorrelated systematic biases in the number of signal and background events in 
each event sample and each detector. Systematic uncertainties are assumed to be completely 
uncorrelated between the neutrino and anti-neutrino runs of the experiment. 

We use the following performance indicators to estimate the sensitivity of the experiment to 
standard oscillation physics 

• 0i3 discovery reach. For each combination of true $13 and true <5cp ; we compute the 
expected experimental event rates, and then perform a x 2 fit assuming a test value of $13 = 0. 
If, for a particular combination of 91™° and 5qp s , the fit disagrees with the simulated data at 
a given confidence level, we say that that this 0^ c and S^p e are within the discovery reach 
of the experiment at that confidence level. In the fit, we marginalize over all oscillation 
parameters except $13 (which is kept fixed at zero) as well as the matter density. 

• Discovery reach for the normal mass hierarchy (NH) For each point in the 9^ e -5^p e 
plane, we simulate the event rates assuming a normal mass hierarchy, and then attempt a 
fit to the simulated data assuming the inverted hierarchy. If the fit is incompatible with the 
data at a given confidence level, we say that the chosen combination of 0\™ c and 5^p e ^ s 
within the NH discovery reach of the experiment. 



• CP violation (CPV) discovery reach. For each point in the 0^ 3 ue -<5Qp e plane, we simulate 
the expected event rates and then attempt fits assuming 5qp = and 5cp = If the fits 
are able to exclude the CP conserving solutions at a given confidence level, we say that the 
chosen combination of &\% c and ^cp e is within the CPV discovery reach of the experiment. 



9 



WC, 200 kt fiducial LAr, 33.4 kt fiducial 




FIG. 3. Sensitivity to standard oscillation physics in a wide band beam for a 200 kt (fiducial) water 
Cerenkov detector (left) and a 34 kt (fiducial) liquid argon detector (right). The results are shown at 3<7 
confidence level. 



• Sensitivity to the octant of 023- For each point in the #23 Ue ~^i3 Ue P iane > we simulate the 
expected event rates and then attempt a fit in which all parameters are marginalized over, 
but 823 is forced to lie in the "wrong" octant, i.e. between 7r/4 and n/2 — #23 UC - K the fit is 
incompatible with the simulated data at a given confidence level, we say that the experiment 
is sensitive to the octant of 823 at that confidence level. 

When discussing non-standard neutrino interactions, we will use the NSI discovery reach as 
a performance indicator, which we define in analogy to the #13 discovery reach: For each set of true 
NSI parameters, we check whether a standard oscillation fit neglecting NSI is compatible with the 
data at a given confidence level. If this is not the case, the chosen NSI parameters are within the 
experimental discovery reach. 



IV. RESULTS 



Standard oscillation 



First, we summarize the physics performance of the default setup, as defined in section III A, with 
respect to three flavor oscillation. Figure [3] depicts the discovery reaches for CPV, #13 and the 
mass hierarchy. In the left hand panel we show results for a 200 kt water Cerenkov detector (WC), 
whereas in the right hand panel we show the corresponding results for a 34 kt liquid argon detector 
(LAr). We have checked that the difference in performance between the 60GeV and 120 GeV 
proton beams, at equivalent power, is very small, and therefore we only show the result for the 
120 GeV beam. It is apparent from this figure that the performance of the two detectors, despite 
a factor of 6 difference in fiducial masses, is quite similar. If sin 2 29\3 < 0.04 the beam upgrade 
provided by Project X, whose results are shown as dotted lines, is a necessity to ensure a mass 
hierarchy determination and to maintain a better than 50% coverage for CP violation. Even for the 
largest possible values of #13 the CP sensitivity would greatly benefit from a luminosity upgrade, 
as also can be seen from figure |j| We have also evaluated the relative precision on sin 2 29\3, defined 
by (sin 2 29f 3 ax — sin 2 28^ n ) / sin 2 26**3 Ue , where 6*™ 11 and Off 1 * denote the lower and upper bounds on 
613 that can be expected for a particular 9\ T 3 ne . We find for sin 2 26*^° = 0.1 that the WC detector 
measures sin 2 29\3 with a relative 3a error between 33% and 39%, depending on the true value of 
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FIG. 4. Ability to determine the octant of 6*23 at 3 a and 90% confidence level for 200 kt WC and and 
34 kt LAr detectors as labeled in the legend. The result is shown as a function of 0| r 3 ue and sin 2 26\ r 3 uc . The 
width of each region is due to the unknown CP phase. 



Sep, while the LAr detector can achieve a precision between 36% and 42%. This result should be 
compared with the accuracy obtainable from reactor neutrino experiments like Daya Bay, which 
will provide a relative error of 18% at 3a C.L. [9]. In figure [4] we study the ability to determine the 
octant of 023 in a WC or LAr detector. The width of the colored regions in the plot is due to the 
marginalization over the unknown CP phase ^p 6, ^or this measurement, we see again that the 
differences in performance between a WC detector and a six times smaller LAr detector are small. 
Determining the octant of #23 is a difficult measurement for any experiment. In particular, for 623 
close to 45° this measurement is only possible for large #13. The asymmetry in sensitivity between 
#23 Ue above or below 45° is due to the partial cancellation between the octant sensitive terms in 
P^e and matter effects. Therefore, the sign of the asymmetry will change if we were to change the 
assumed true mass hierarchy from normal to inverted in our calculation. 

B. 2 nd maximum 

Now that we have established the baseline performance, we can turn our attention to the central 
question of this paper: what is the quantitative impact of data from the 2 nd oscillation maximum 
on the physics sensitivities? This question specifically neglects the issue of how the robustness 
of an experiment with respect to unforeseen systematical effects improves due to the data from 
2 nd oscillation maximum. However, the current analysis does include known systematic effects 
like normalization errors of backgrounds and signal. Obviously, if the data from the 2 nd oscillation 
maximum can be collected at no or only very small cost, we are well advised in using it, even if only 
to check whether our assumptions about the performance of the experiment and the underlying 
physics model are correct. However, in case that obtaining this data turns out to be costly, we 
need to understand in a quantitative way how much one would lose by not having it. 

For the baseline setup discussed in the previous section, we can show that for all standard 
oscillation measurements, with the exception of the sensitivity to the mass hierarchy, there is 
virtually no difference between an analysis which includes both maxima and one where we ignore 
all data with energies below 1.25 GeV. For the mass hierarchy measurement, the improvement 
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as labeled in the legend of each panel, on the scaling 



Shown are lines of constant CP fraction in the X2~sin 26**™° resp. 



X2-|(e^? 1 ) e | plane. The results are given at 3 a confidence level. 



happens for those values of the CP phase where the 7r-transit phenomenon (see section II A ) would 
strongly reduce the sensitivity. Even there, the data from the 2 nd maximum is statistically not 
significant enough to improve the sensitivity to the level it would have if there were no 7r-transit. 

This result does not imply that the 2 nd maximum makes no quantitative difference at all, it just 
shows that, with the specific experimental setup chosen, the data sample in the 2 nd maximum is 
too small and the backgrounds are too large (see figure [2J in order for that data to make a sizable 
contribution to the overall x 2 - Therefore, we will now study how the sensitivities change if there 
are more events in the energy region below 1.25 GeV. If we just were to scale up the number of 
events in the 2 nd maximum, obviously, we always would find that the x 2 becomes larger, since 
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it is a monotonically increasing function of the total number of events. However, at constant 
beam power, the number of pions produced in the target is constant as well. Therefore, any 
beam optimization is equivalent to selecting a different subset of pions leading to different neutrino 
spectra. Assuming, furthermore, that the acceptance of the horn and beam pipe are finite and 
fixed, any optimization is just a reshuffling of pions and hence neutrinos of different energies, with 
the total number of neutrinos remaining fixed. This inspires the following parametrization: Let </> 
be the total flux of Vn in the beam, and let (j>\ (^2) be the partial flux in the energy window above 
(below) 1.25 GeV, corresponding to the 1 st (2 nd ) oscillation maximum. To assess the importance 
of the second maximum to the experimental sensitivity, we vary the fraction of neutrinos below 
1.25 GeV, while keeping the total flux constant. More specifically, we scale 02 with an efficiency 
factor X2, and <\>\ with an efficiency factor x\ = [cpi + (1 — X2)4>2]/<j>i, so that 0i + 02 does not 
change. X2 can take values between and 0/02, with X2 = 1 corresponding to the setup defined 



in section III A For each fixed X2, we compute x\ separately for the neutrino beam and the anti- 



neutrino beam. To summarize this parametrization, we have 



+ 4>2 and 4>2 — > X2 4>2 and <j>\ — > — %2 02 and X2 G 



0. 



(10) 



In the first three panels of figure [5| we show the sensitivity to standard oscillation physics as a 
function of X2- We see that for all performance indicators except the mass hierarchy, the optimum 
occurs at X2 < 1, which implies that we rather have more events in the 1 st maximum instead of 
sharing them with the 2 nd maximum. 

An altogether different way to access events in the 2 nd oscillation maximum is to use a second 
detector at a longer baseline. If the two detectors are to be in the same beam and the first one is 
on-axis, the second one has to be off-axis due to the curvature of the Earth's surface. Alternatively, 
one can also imagine scenarios in which both detectors are off-axis, either with identical off-axis 
angles or with different ones. For example, in the T2KK setup |38l I45j . it has been proposed to use 
Super-Kamiokande (or a larger Water Cerenkov detector at the same site) as the detector sensitive 
to the 1 maximum, and supplement it with a second Water Cerenkov detector at a baseline of 
about 1,000 km on the east coast of Korea. Another proposal [36] puts a liquid Argon detector at 
about 600 km on the island of Okinoshima. In both cases, due to the different off-axis angles, the 
second detector will be predominantly sensitive to the 2 nd maximum. A superficial comparison of 
the obtainable sensitivities indicates a similar physics performance, where most of the differences 
is attributable to the different overall exposure [42] . In order to allow for a direct comparison with 
the results derived in this paper, we refrain from comparing these setups in detail and study instead 
the effects of the addition of a second baseline to the setup we have introduced in section |III A 



Since we are interested in the question of the general impact the 2 nd maximum can have, we will 
neglect the actual geometry and assume that we have two identical beams, which allows us to 
put the second detector on-axis into this second beam. This is clearly an unrealistic and overly 
optimistic assumption. It amounts to doubling the number of protons on target and in contrast 
to the proposals centered around Super-Kamiokande, which all exploit a single beam, leads to 
higher event rates in the 2 nd maximum due to it being accessed in an on-axis beam. However, it 
allows us to estimate the maximum effect that events from the 2 nd maximum could have under ideal 
circumstances. In other words, if we do not observe an overwhelming increase in performance under 
these most favorable conditions, then we can safely conclude that the 2 nd oscillation maximum, 
despite its theoretical merits, in practice is not useful in a superbeam experiment. The results of 
this analysis are shown in figure [6j For the measurement of #13 and CP violation, the performance 
optimum occurs for a detector location very close to 1 , 300 km, which is the position of the first 
detector. The sensitivity to the mass hierarchy shows a strong preference of baselines around 



2, 500 km, but we remark that a similar effect is also seen with only one detector, see figure 10 and 
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also reference [39]. Interestingly, at this distance the first oscillation minimum is at the peak of 
the beam flux and thus both the 1 st and 2 nd maximum contribute about equally to the rate. 

A further question is whether information from the 2 nd maximum can help in determining the 
octant of #23- To this end we computed the sensitivity to the octant in the same way as shown in 
figure 4 but constraining the data to the 1 st maximum only. The results are identical to the one in 
figure 4; thus, we find that data from the 2 nd maximum does not improve the senstivity to discern 
the octant of #23- 

The conclusion for three flavor oscillation in this case is the same as with only one detector: 
The 2 nd maximum does not help with the measurement of #13 or the CP phase, but it enhances the 
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FIG. 7. Discovery reach for neutral-current-like NSI in a WC detector (colored bars and arrows), compared 
to current model-independent limits on the different parameters [3, 28 , 29J (gray shaded regions) . The length 
of the colored bars is due to the unknown phase associated with the new interaction. 

ability to measure the mass hierarchy. The gains in mass hierarchy sensitivity could be substantial 
under favorable conditions, but are only moderate in practice. 

C. Non-standard Interactions 

Next we turn our attention to the question whether the 2 nd maximum is useful for new physics 
searches. For a given set of oscillation parameters the relative strength of the signals in the 1 st 
and 2 nd maximum is well understood within the standard three flavor oscillation framework, and 
therefore any deviation should stem from new physics. In order to be able to perform a quantitative 
analysis of this problem, we will restrict the new physics to the form of non-standard interactions, 



their underlying physics and parametrization have been described in section II B 

In figure [7] the discovery reach for neutral current like NSI is shown for our standard setup. For 
each set of bars, only one NSI parameter was allowed to be nonzero at a time, i.e. we do not include 
correlations between different NSI parameters. The length of the bars is due to the unknown phase 
of the non-standard parameters, whereas the different colors are for different subsets of the data as 
explained in the legend. The gray shaded areas indicate the current model independent bounds on 
these parameters |3 ESI 12210 

in cases where there is no gray shaded area, the current bounds are 
of order one. Note that possible correlations between and other parameters are equivalent to 



We have converted the 90% C.L. bounds given in refs. {51 l28j 129] to the 3cr confidence level assuming Gaussian 
errors. 
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correlations with the matter density, which are included in our simulations by the matter density 
uncertainty. Also, correlations between and the [i-t sector (e^, e™ T , e™ ) are small because 

e™ and affect mainly the appearance channel, while e^, e™ T and are most relevant in the 
disappearance channel, see e.g. references |26} 127]; this has been shown for the £™-£™ correlation 
by explicit numerical calculation in Ref. [47J. The strongest improvement in bounds happens for 
flavor- changing NSI, but this improvement is hardly dependent on the data from the 2 nd oscillation 
maximum. The results for a liquid argon detector are very similar and lead to the same conclusion. 

Based on our somewhat negative result with respect to the 2 nd maximum in the case of standard 
three-flavor oscillation, we can study the effect of an artificial enhancement of statistics in the 
2 nd maximum also in the presence of non-standard interactions. The result of rescaling the flux 
according to equation 10 is shown in the lower right hand panel of figure [5] for the case of NSI in 



the e— \i sector. The performance optimum occurs at X2 < 1, which implies that the 2 nd maximum 
is not useful in this case. The option of using a second detector to access the 2 nd maximum exists 
also for NSI studies and the result is shown in the lower right hand panel of figure |6] With a second 
detector between 3, 000 — 6, 000 km the sensitivity would be improved by less than a factor of two. 
The wide baseline range over which this improvement happens makes it seem unlikely that this is 
entirely due to the 2 nd maximum. Also, in contrast to the case of the mass hierarchy measurement, 
where there was an improvement both for rescaling the flux and considering a second detector, here 
we see improvement only for the second detector, which points to overall increased matter effects, 
standard and non-standard, as the source of this improvement. Qualitatively a similar improvement 



is obtained by just using one detector at a longer baseline as shown in figure 10 In any case, the 
improvement is relatively moderate. 

The previous statement about the relative unimportance^] of correlations between NSI parame- 
ters does not hold if one allows for the simultaneous presence of non-standard effects in the neutrino 
production, propagation, and detection processes. In particular, in this case the so called confusion 
theorem obtains |48| : Assume that #13 = 0, but there are charged current NSI between v T and 
electrons in the detector [e d re ^ 0), and neutral current NSI between v e and v T in the propagation 
(e™ / 0)0 If we furthermore assume that the parameters obey the relation 

4e = re? T , (11) 

with r being an order one parameter determined by whether the NSI couple to quarks or leptons, 
then the event rate spectra for both neutrinos and anti-neutrinos in the — > v e channel are the 
same as for standard three-flavor oscillations with 

sin^ 13 = r 2 (C) 2 l + C0 2 S2g23 . (12) 

This is the confusion theorem. Subsequently, it was discovered that for sufficiently high beam 
energies, muons from the decay of r from v T charged current interactions can be used to resolve 
the confusion at least for parts of the parameters space [49]. For the setup considered here, the 
average beam energy is close to m T , and therefore T-production in charged current interactions 
from v T will be strongly suppressed, and therefore no muons from r-decays will be observed. Thus, 



the confusion theorem should apply. Still, it is important to note that equations 11 and |12| were 
obtained from a perturbative expansion of the oscillation probability. This expansion is strictly 
valid only for energies around the 1 st maximum. Thus the question arises to which degree the 
confusion theorem applies to the wide band scenario considered here. 



3 While there may be correlations between the various parameters, the fact that all relevant e involving /i-type 
flavor are tightly constrained, should make the correlations practically negligible. 

4 The original confusion theorem was derived in the context of a neutrino factory, where the appearance signal stems 
from v e —¥ oscillations. There, it is nonzero together with a CC-like NSI in the source (e^ T 7^ 0) which 
causes the confusion. Here, we are considering the T-conjugate oscillation channel, and hence we need a CC-like 
NSI in the detection process instead. 
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corresponding results if one assumes #13 to be constrained by a measurement at Daya Bay. 



A partial answer is shown in figure |8j where the sensitivity to sin 2#i3 , CP violation and the 
mass hierarchy is depicted for various levels of NSI. The different lines are obtained by allowing 
successively larger values of |e™| and |e^ e | in the fit, as indicated in the plots. The largest values 
used in figure [8] correspond to the current 3 a bounds according to references |28|, [29] j^] In the left 
panel (#13 discovery reach) and in the right panel (discovery reach for the normal mass hierarchy), 
we allow and e^ e to be complex with arbitrary phases, while in the middle panel (discovery 
reach for CP violation) , we assume the phases of the NSI parameters to be or ir since we want 
to consider only CP conserving solutions in the fit. The rightmost line in the left hand panel of 
figure [8] confirms the validity of the confusion theorem in equations 11 and 12 We indeed observe 
a deterioration of the sensitivity by nearly an order of magnitude in sin 2 26\™ c . However, at the 
same time we see that alone accounts for most of this deterioration since the difference between 
the rightmost line and the line next to it is relatively small. Thus, the confusion theorem seems 
to apply in essence, but in practice the small number of events around the sensitivity limit does 
not require the presence of NSI both in propagation and detection, because spectral information 
is not statistically significant. For the same reason, the information from the 2 nd maximum plays 
no role, since at the sensitivity limit the event sample from the 2 nd maximum is statistically not 
significant. We have performed the same scaling analysis as presented in figure [5] also in this case 
and find that the 2 nd maximum is not useful in controlling the effects of NSI in propagation and 
detection. 

In the middle panel of figure [8] we show the impact of NSI on the ability to discover CP 
violation, and the result exhibits the same qualitative features as the one for the discovery of #13. 
At a quantitative level, this measurement is more sensitive to the deleterious effects of NSI since 
it relies on smaller, more difficult signatures also in the standard oscillation case. Therefore, we 
observe a complete loss of discovery potential for values of the NSI an order of magnitude below 
the current bounds. One may speculate that using a precision measurement of $13 by the Daya 
Bay experiment could mitigate the correlation with NSI, however, as the dashed lines conclusively 
demonstrate, this is note the case. Thus, at the current level of analysis we are forced to conclude 



Again, we have converted 90% C.L. limits to 3 a constraints assuming Gaussian errors. 
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that the ability to discover leptonic CP violation in the next generation of superbeam experiments 
is not robust with respect to the presence of new physics. Neither a #13 measurement by Daya 
Bay nor events from the 2 nd maximum can resolve this problem. A direct measurement of in 
a short-baseline neutral current neutrino scattering experiment to improve the upper limits is out 
of question, since this would require to measure the flavor of the outgoing neutrino. As shown in 
figure 4 of reference |48j , combining different baselines does not affect the validity of the confusion 
theorem and therefore does not present a viable strategy to address the confusion problem either. 

In the rightmost panel of figure [8] we show how much the discovery reach for the mass hierarchy 
is diminished, and the result is equally drastic as in the case of CP violation. 

A near detector can help to improve bounds on NSI in the neutrino production and detection 
processes, but our results show that NSI in propagation alone are sufficient to cause serious prob- 
lems for the measurement of the standard oscillation parameters. Besides, a direct measurement 
of e^ e is a difficult proposition for a number of reasons. First, there are no v T in the beam, so the 
only way of constraining e^ e would be to constrain e s eT instead and to make use of the fact that the 
two parameters are usually related since they can both arise from the same non-standard coupling 
between two light quarks, an electron, and a v T [26J. However, this relation between e^ e and e s eT 
is not model-independent. For instance, a parity-conserving non-standard operator can lead to 
nonzero e!~ e , but will not contribute to neutrino production in pion decay, so that e s eT = [50 J . 
Moreover, even measuring e s eT is very challenging because the initial flux of electron neutrinos in 
an LBNE-like beam is very small, less than 1%, the kinematic suppression of r-production is large 
with the available beam energy, and r-identification is notoriously difficult and typically has a low 
efficiency. For these reasons, we conclude that near detectors will not solve the problem of possible 
confusion between standard oscillations and NSI. The specific setup considered here is in some 
sense a best case scenario, since it has relatively high statistics, a lot of spectral information, and 
makes use of two oscillation maxima. There is no reason to expect that experiments like T2K and 
NOfA will be less affected by the confusion problem, quite the contrary, as shown in reference j26j . 



V. SUMMARY AND CONCLUSIONS 



The goal of this paper is to quantitatively understand the benefits, or lack thereof, of studying 
two oscillation maxima simultaneously in a long-baseline neutrino oscillation experiment. To this 
end we have chosen a specific example of experimental setup, which closely resembles the current 
plans for the Fermilab-DUSEL Long Baseline Neutrino Experiment (LBNE). There are two ways 
to access the 2 nd oscillation maximum: either using a broad neutrino energy spectrum to cover the 
1 st and 2 nd maximum in the same detector, or using two detectors at different baselines in the same 
beam. In LBNE, the natural method is to use the same detector and a wide energy spectrum. For 
this approach we found that there is no apparent benefit from using the 2 nd maximum (figure [5]). 
This remains true even if it were possible to shift a larger portion of the total neutrino flux into 
the energy range of the 2 nd oscillation maximum. The measurement of the mass hierarchy does 
improve with events from the 2 nd maximum, but the improvement is limited, and since it would 
come at the cost of losing events in the 1 st maximum, which in turn negatively impacts the other 
measurements, a trade-off between these conflicting requirements has to be found. We have also 
investigated the option of using the same beam but two detectors at different baselines, which 
is the natural option for extensions of T2K [38, 45, 46J. The results are similar to the previous 
case for the measurement of #13 and CP violation, while the improvement in the sensitivity to the 
mass hierarchy is somewhat stronger in this case due to larger matter effects at the longer baseline 
(figure [6]) . As far as the possible detection of new physics — parametrized here in the framework 
of neutral current non-standard interactions (NSI) — is concerned, the sensitivity improves slightly 
for a second detector at a longer baseline. Note that neutral current NSI measurements prefer 
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longer baselines in general since they essentially correspond to new neutrino matter effects, and 
matter effects are larger at long baseline. Therefore a similar improvement can also be observed 
for a single detector setup with a longer baseline (figure 10). 



Finally, we have revisited the so called confusion theorem, which has been discovered in the 
context of neutrino factories |48j . The confusion theorem states that certain combinations of 
charged and neutral current NSI modify the neutrino oscillation probability in the same way as a 
nonzero value of #13 does. In particular, the effects of the NSI parameter e™ are problematic since 
this parameter is only weakly constrained by an 0(1) bound. In other words, in the absence of 
sufficiently strong bounds on NSI, it is very hard to establish a bound on #13. In the context of 
neutrino factories, muons from r decays have proven to be a loophole in the confusion theorem |4i)J , 
but in a superbeam experiment the beam energy is so low that T-production is kinematically 
suppressed and the arguments from reference |49j do not apply, as has been shown in the context 
of T2K and NOi^A [26]. In the present work we have extended those earlier results to superbeam 
experiments which have events from the 1 st and 2 nd oscillation maximum. We have found that 
the confusion theorem holds (figure [8| and the sensitivity to sin 2 28\^ deteriorates by one order 
of magnitude if the possibility of NSI is taken into account. Data from the 2 nd maximum has no 
effect on this conclusion since it is statistically not significant enough. Moreover, we have shown 
that even the effect of neutral current NSI (e™ 7^ 0) alone is sufficient to seriously impact the 
sensitivity to sin 2 2^13. This is an especially problematic limitation since cannot be constrained 
by a near detector measurement. For CP violation measurements, the possibility of complex e™ 
has been considered and the impact is dramatic: Even for an order of magnitude below the 
current bound a complete loss of sensitivity ensues. Reactor experiments will not be affected by 
e™ and thus can provide a clean measurement of sin 2 2#i3 , but we have shown that even precise 
knowledge of #13 from a reactor experiment does not solve the problem. The presence of data from 
the 2 nd maximum does not provide immunity from the confusion theorem and it has an overall 
very small quantitative impact. This remains true even if the majority of the neutrino flux were 
shifted into the 2 nd maximum. In comparison, T2K and NOi^A have smaller statistics and observe 
a much smaller range in L/E. Therefore, the impact of the confusion theorem is more severe in 
these experiments [26]. 

To put our conclusion on the possible impact of large NSI into perspective, we should emphasize 
that NSI large enough to be problematic for T2K, NOz/A, or LBNE, are not a generic prediction 
of extensions of the Standard Model. In particular, from a model-builder's point of view, new 
particles at or above the electroweak scale are rather unlikely to have a sizable effect on the 
neutrino sector [244 I25j . On the other hand, new physics at a low scale may still lead to large 
NSI [221 [23], and since neutrino physics has taught us in the past that theorists' prejudices may 
be wrong, one cannot discard such possibilities. 

In summary, we find that only the determination of the mass hierarchy benefits slightly from 
using the 2 nd oscillation maximum in a long baseline neutrino oscillation experiment. All other 
measurements, including non-standard neutrino interactions, are not improved compared to the 
case where only events from the 1 st maximum are used. We confirm that the "confusion theorem" , 
which states that certain types of NSI can mimic the effects of nonzero #13, remains valid even 
if data from the 2 nd maximum is available. In particular, we have shown that the presence of a 
non-standard coupling between electron neutrinos and r neutrinos with a complex coefficient 
can completely destroy the sensitivity to CP violation and the mass hierarchy. Since it is known 
from the literature [H] that combining different baselines does not alleviate this problem, and we 
have explicitly shown that adding reactor neutrino data does not work to this end either, it seems 
that there are no simple remedies for the confusion problem. 
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Appendix A: Possible alternative setups 



In this appendix we study simple variations of the basic experimental setup considered in this 
paper. The results given here will allow to extrapolate the effects of changes in the exposure 
and baseline. In figure [9] we show the discovery reach for standard oscillation as well as NSI as 
a function of exposure. Obviously, the higher the exposure the better the sensitivity. Discovery 
of CP violation has the largest demand for high exposure and conversely is most at risk if the 
luminosity were to turn out smaller than expected. 

Figure 10 illustrates the dependence of the discovery reaches for #13, mass hierarchy, CPV 
discovery and the NSI parameter as a function of the baseline. The lines and corresponding 
shades are iso-contours of CP fraction. The optimum occurs for CP violation and the discovery 
of 0i3 around 1,500 km, whereas the optimum for the discovery of the mass hierarchy and is 
around 2, 200 km. This result confirms that L = 1, 300 km is a reasonable comprise for the neutrino 
beam assumed in this paper, and slightly longer baselines around 1,600 km would perform very 
similarly. 
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